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Introduction to hydrodynamic 
modelling and granular dynamics 

The present text introduces hydrodynamic modelling principles in the context of batch wet granulation 
and coating systems and it reviews the latest achievements and proposals in the scientific literature in this 
field. The text concerns primarily the Eulerian and the Lagrangian modelling technique. In accordance 
with some of the latest published Ph.d. thesis in the field of hydrodynamics modelling, the Lagrangian 
technique is divided into a soft-particle and a hard-sphere approach. The text further presents some 
of the latest trends and results from the growing field of applying Computational Fluid Dynamics and 
Discrete Element Modelling in the field of modelling fluid bed granulation processes. Further, a number 
of other granule dynamic modelling principles including the Finite Element and Monte Carlo techniques 
are introduced. 

The text is aimed at undergraduate university or engineering-school students working in the field 
of mathematical or chemical and biochemical engineering. Newly graduated as well as experienced 
engineers may also find relevant new information as emphasis is put on the newest scientific discoveries 
and proposals presented in recent years of scientific publications. It is the hope that the present 
introductory text will be helpful to the reader - particularly in the early stages of the process of working 
with hydrodynamics in a granulation context. The comprehensive literature list may also hopefully be 
an inspiration for further reading. 

I alone am responsible for any misprints or errors and I will be grateful to receive any critics and/or 
suggestions for further improvements. 


Copenhagen, September 2006 
Peter Dybdahl Hede 
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1 Hydrodynamic modelling and 
granular dynamics in respect to 
fluid bed processing 

As it has been emphasised elsewhere (e.g. Hede, 2006b), a lot of research in fluid bed processing 
has focussed on modelling and understanding of the separate growth mechanisms associated with 
agglomeration. Much of this understanding of these separate mechanisms has been integrated into 
population balance models as it was presented in Hede (2006c). The available population balance models 
almost never take into account the hydrodynamic properties and influence on the fluid bed process, 
although the nature and detailed knowledge of fluidisation must be regarded as a prerequisite for precise 
modelling - especially when it comes to scaling of fluid bed processes. Hydrodynamic modelling of 
fluid bed systems seeks to include the effect of fluidisation although this approach is somewhat different 
from the population balance approach 1 . A hydrodynamic approach to fluid bed systems is a fairly new 
discipline that has developed in the last five to ten years due to the enormous increase in computer 
processing power and algorithm development. It is nevertheless a difficult but promising discipline 
being a prerequisite for discrete element modelling and computational fluid dynamic modelling which 
are both subjects to be presented in the present chapter. 

Hydrodynamics is fluid dynamics applied to liquids, and as the bed load is fluidised in the fluid bed 
vessel during processing, hydrodynamics is often used as a term rather than the more general term fluid 
dynamics 2 . As with other typical fluid dynamic problems, a fluid bed hydrodynamic problem often 
involves the calculation of various properties for the fluidised particles such as velocity, pressure, density 
and temperature as function of space and time. The system consisting of the fluidised solid particles 
and the fluidising gas is often treated as a two-phase gas-solid flow system (Goldschmidt, 2001 and 
Goldschmidt et al., 2003). Modern approaches by Goldschmidt (2001) reduce this simplification by 
treating the fluid bed system as a multiphase system, which allows a detailed hydrodynamic modelling of 
the fluid phase of top-spray fluid beds. The description of such systems seeks in any case to account for 
the inherent complexity of dense gas-particle flows, which in its turn can be related to particle-particle 
and particle-wall interactions as well as gas-particle interactions (Goldschmidt, 2001). Hydrodynamic 
modelling is rarely combined with mass transfer or chemical/physical reactions and mechanisms as first 
attempts by Samuelsberg & Hjertager (1996) strongly indicate the prior need for valid and well-proven 
hydrodynamic models. 
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On an overall scale, there are basically two types of hydrodynamic models being the Eulerian models 3 
and the Lagrangian models of which the first is commonly known as Computational Fluid Dynamic 
models (CFD) and the latter is commonly referred to as Discrete Element Models (DEM) 4 . Both types 
consider the gas phase as a continuum but there are quite a few differences in the modelling approaches 
and assumptions. The following sections will briefly introduce the basic principles and differences. An 
in-depth treatment of fluidised bed hydrodynamics and granular dynamics of two-phase flows should 
be found in Hoomans (1999) and Goldschmidt (2001). 

1.1 Eulerian models 

In Eulerian models the gas and the solid phases are treated as interpenetrating phases, and the theory 
behind such models is basically an extension of the classical kinetic theory that takes non-ideal particle- 
article collisions and gas-particle drag into account (Goldschmidt, 2001). In this scheme, collections of 
particles are modelled using continuous medium mechanics. The solid particles are generally considered 
to be identical having a representative diameter and density, meaning that the particle phase is volume 
averaged (Taghipour et al., 2005 and Depypere, 2005). The general idea in formulating such a multi¬ 
fluid model is to treat each phase as an interpenetrating continuum and therefore to construct integral 
balances of continuity, momentum and energy for both phases with appropriate boundary conditions 
and jump conditions for the phase interfaces. Since such a resulting continuum approximation for the 
solid phase has no equation of state and obviously lacks variables such as viscosity and normal stress, 
certain averaging techniques and assumptions are required to obtain a momentum balance for the solid 
phase (Pain et al., 2001). 

Although constitutive relations according to the kinetic theory of particle flow have been incorporated 
into recent models (e.g. Gidaspow et al., 2004 and Chiesa et al., 2005), pure CFD models for fluid bed 
granulation still suffer from the fact that the contact between fluid, particles and boundary surfaces is not 
considered explicitly with respect to particle inertia and the mechanical properties of the particles. This 
limits the ability of CFD multiphase models to adequately represent particle-particle and fluid-particle 
interactions thereby reducing the accuracy of the prediction of both the fluid and the particle dynamics 
(Fan et al., 2003 and Sun, 2002). This inaccuracy can be overcome by explicit calculations of the particle 
contact mechanics in a particle-scale reference frame using a Lagrangian approach as it will be presented 
below. Considering the required computational power and complexity, gas-particle flow fields calculated 
with the multi-fluid interpenetrating approach of the Eulerian granular multi-phase model is still a fast 
method to calculate flow fields, as it is well known from simple particle systems as spray-drying and 
conveying systems etc. Due to the obvious need for accounting precise particle level properties into fluid 
bed hydrodynamic models, pure Eulerian CFD models 5 must be regarded as inappropriate even in an 
industrial context. Hydrodynamic fluid bed scaling attempts on empty vessel basis are often based on 
Eulerian models (e.g. Krishna & van Baten, 2001 and Cooper & Coronella, 2005) but the models often 
turn out to be inaccurate when particles are being processes. Recent attempts by e.g. Depypere (2005) 
and Lettieri et al. (2003) quite clearly illustrates that future development in the field of hydrodynamic 
fluid bed modelling should be concerned with Lagrangian principles rather than Eulerian. 
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1.2 Lagrangian models - Discrete Element Methods 

The Lagrangian approach maybe seen as an extension of the Eulerian, in that Lagrangian model describes 
the solids phase at particle level and the gas phase as a continuum. In the two-phase flow situation the 
Newtonian equations of motion for each individual particle are solved with inclusion of the effects of 
particle collisions and forces acting on the particles by the gas (Goldschmidt, 2001). That is, Newtons 
law of motion is simultaneously solved for a large number of particles either in a computational unit 
cell with periodic boundaries or on a computational domain representing the entire fluid bed vessel or 
its subset. In general the following equations are solved (Cameron et al., 2005): 


dV i 17 

m, — L = F 
1 dt 


and 


( 1 . 1 ) 


j d®i 

1 dt 


= M; 


( 1 . 2 ) 


in which m. is the particle mass, v. is the velocity vector, t is time, I. is the moment of inertia, co. is the 
angular velocity vector, M. is the net torque vector and F. is the net force vector acting on particle i. The 
net force vector F. is often written as the sum of three contributions (Cameron et al., 2005): 


F = F h + F p + F e 


(1.3) 


where F H is the force due to fluid-particle interactions also known as the drag force, F p is the force due to 
particle-particle interactions during collisions and F E is a force acting on the particle due to an external 
field being e.g. the gravitational field. F E is sometimes completely ignored and F H is often estimated 
from empirical equations, and many discrete element approaches involves different assumptions and 
expressions for the net force vector in equation 1.3 (Goldschmidt et al., 2004). 


Besides obviously being far more precise than the Eulerian models, such discrete particle models do 
not require additional closure equations for the suspended particulate phase since they compute the 
motion of every individual particle, taking collisions and external forces acting on the particles directly 
into account (Goldschmidt, 2001). The Lagrangian approach may roughly be divided into two groups 
being the soft particle and the hard-sphere approach, both of which will be briefly introduced below. 
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1.2.1 Hard-particle models 

In hard-sphere models the particles are assumed to interact through instantaneous, binary collisions. 
A sequence of collisions is processed one collision at the time in order of occurrence (Hoomans, 1999, 
Hoomans et al., 2000 and Tsjui et al., 1993). Hard-sphere models are also referred to as event driven 
models since a sequence of collisions is processes in which all particles are moved until the next collision 
occurs. Particle collision dynamics are described by collision lays, which account for energy dissipation 
due to non-ideal particle interaction by means of the empirical coefficients of normal and tangential 
restitution and coefficient of friction (Goldschmidt, 2001). The dissipative particle interaction in particle 
media makes these systems significantly different from molecular systems where energy associated with 
collision always is conserved. This means that energy has to be continuously supplied to the particle 
system in order to keep the particles in motion. This can for instance be achieved by applying a shear 
rate through proper choice of boundary conditions as suggested by Campbell & Brennen (1985). It has 
become common to choose the collision particle partners and sequences based on the relative approach 
velocity (Hoomans et al., 1996 & 2000). 

Based on work by Hoomans (1999), Goldschmidt (2001) and Goldschmidt et al. (2003) developed a 
hard-sphere discrete particle model for gas-fluidised beds which captures the principles of basic two 
dimensional hard-sphere modelling well. This model computes the motion of every individual particle 
as well as droplet in the system considering the gas phase as a continuum. Micro-scale processes such 
as particle-particle collisions, droplet-particle coalescence and agglomeration are taken into account by 
simple closure models. Distinction is made between three types of entities being dry particles, wetted 
particles and droplets. All three types are assumed to be spherical, and encounters are detected as soon 
as contact occurs at a point on the line joining the centres of the two entities. In addition, six types 
of encounters are distinguished: encounters among dry particles described by hard-sphere collision 
laws from mechanics; Droplet-droplet encounters described by hard-sphere collisions laws as well, as 
they are assumed to be repulsive for atomised liquid droplets with a typical radius of 50 pm, colliding 
with small mutual differences; Encounters between droplets and dry or wetted particles, described as 
coalescence; Encounters between dry or wetted particles and a wall, described by hard-sphere collision 
laws; Encounters between droplet and walls, resulting in removal of the droplet from the simulation and 
last; Encounters between a wetted particle and another particle, leading to either rebound described by 
hard-sphere collision laws or agglomeration. Which of the two situations that occurs depends on the 
odds of the particles hitting each other on a wet spot (Goldschmidt et al., 2003). It is further assumed 
that a new particle entity is formed at the position of the centre of mass of the original entities upon 
coalescence as it is sketched in figure 1. 
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Figure 1: Repositioning and merging of particles upon coalescence or 
agglomeration (Goldschmidt et al., 2003). 

Mass, momentum and volume of the new entities are conserved and transferred to the newly formed 
particle, where after the original entities are excluded from the simulation. In case of coalescence, the 
area on the newly formed particle covered by liquid depends on the original particle size, the size of the 
droplet and a defined minimum liquid layer thickness as sketched in figure 2. 
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Figure 2: Liquid layer formation upon coalescence (Goldschmidt et al., 2003). 


In case of agglomeration, the wetted area available for subsequent agglomeration is reduced by the 
projected area of the smallest particle, to account for liquid bridge formation and the masking of the 
wetted surface, which cannot be reached anymore because the newly agglomerated particle is in the 
way This is sketched in figure 3. 


O 
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Figure 3: Masking of wetted surface for subsequent agglomeration 
(Goldschmidt et at, 2003). 

Although inclusion of liquid and gas inside the pores is taken into account for agglomerates containing 
more than three primary particles in the model by Goldschmidt et al. (2003), the model does however not 
account for particle deformation, liquid spreading, breakage of droplets and agglomerates, and is further 
limited by the two-dimensional geometry Even so, the simulation using the simple two-dimensional 
Goldschmidt model is only possible for 50,000 granules at the time thereby being only comparable 
to experimental data from laboratory-scale fluid beds. This clearly illustrates the need for increased 
computational power if the hard-sphere principles should be used for any industrial purposes, and also 
indicates why the hard-sphere approach has first been used within the last ten to fifteen years although 
the principles were introduced originally in the late fifties. In recent modelling attempts, hard-sphere 
models are mainly concerned with rapid particle flow 6 (Lian et al., 1998) and the majority of the latest 
discrete element attempts concerns the soft-particle approach. 
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1.2.2 Soft-particle models 

The soft-particle approach differs from the hard-particle approach in that it treats interparticle collisions 
as a continuous process that takes place over a finite time. In such models the particles are assumed 
to undergo deformation during their contact, where the contact forces are calculated from a simple 
mechanical analogy involving a spring, a dashpot and a friction slider as the normal and tangential 
component of forces are expressed as the sum of forces due to the springs and dashpots, and the normal 
and tangential velocities are expressed in terms of the relative velocity prior to collision (Gera et al., 
1998). Such inter-particle bond models are particular suitable for the modelling of impact breakage of 
pre-existing agglomerates which undergo some sort of brittle fracture (Thornton & Cismocos, 1999). 
This is why the majority of discrete element simulations of agglomerate strength use the soft-particle 
approach as mentioned in the previous section. The principle of the linear spring-dashpot model is 
sketched in figure 4. 



Figure 4: Contact force model for soft particle modelling (Gera et al., 1998). 

With soft-particle simulations, the interactive forces exerted on each particle are computed as continuous 
functions of the distance between contiguous particles and are based on physically realistic interaction 
laws. Soft-particle models are also referred to as time driven models as all particles are moved over a 
certain period of time where after the collision dynamics are computed from the particle overlaps. In 
case a particle is in contact with several other particles, the resulting contact force follows from the 
addition of binary contributions (Goldschmidt, 2001). Compared to the hard-sphere principle this 
approach is computationally intensive and requires even higher computational demands than the hard- 
sphere simulations, but does as a clear advantage provides information on the structure and dynamics 
of particular materials including details of positions, velocities, forces and energy partitions (Lian et al., 
1998). This makes soft-particle modelling useful in the simulation of the deformation and breakage of 
agglomerates. In soft-particle simulations of coalescence, agglomerates are modelled as assemblies of 
primary particles, which often are assumed to be spherical and elastic. 
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In a representative example of a soft-particle coalescence model, Lian et al. (1998) developed a model 
in which each agglomerate comprised 1000 randomly packed primary particles with the interparticle 
interactions modelled as the combination of the solid-solid contact forces and the principles for the 
pendular liquid bridges presented in Hede (2006a). It is assumed that the binding liquid present at 
particle junctions completely wets the particles. With this set-up, the coalescence of two randomly packed 
agglomerates each consisting of 1000 primary particles was studied at different impact velocities. Initially 
1000 primary particles were randomly generated in each of two specified spherical regions that were not 
able to touch. Centripetal gravity fields were then applied to the two spherical regions in order to bring 
the particles together and, depending on the size of this field, the regions would impact each other at 
different velocities. When two primary particles collided, a pendular liquid bridge was assigned to that 
contact. The kinematic energy was eventually dissipated due to the retarding effect of the viscous liquid 
bridges. When equilibrium was achieved, the centripetal gravity was removed. Pendular liquid bridges 
then held the particles together. The impact simulation was implemented for a range of initial relative 
velocities between the two agglomerates and for different interstitial binder fluid viscosities, and for 
each collision, the simulation was continued until a major proportion of the initial kinetic energy was 
essentially dissipated. Examples of computer simulated wet agglomerates after impact at a velocity of 
2.0 m/s with different interstitial viscosities ranging from 1 mPa s to 100 mPa s can be seen in figure 5. 
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Figure 5: Computer simulations of agglomeration. 

Visualisations of computer simulated pendular state wet agglomerates after impact at a velocity of 2.0 
m/s for interstitial viscosities of: (a) 100 mPa s, (b) 10 mPa s, (c) 1 mPa s (Lian et al., 1998). 

The results by Lian et al. (1998) illustrate the possibilities regarding the simulation of the coalescence 
situation at particle level, but with the large computer processing time even for the binary collision situation 
it is obvious that such an approach cannot easily be extended to model the high number of particles inside 
a commercial fluid bed. Even if it at some point will be possible due to the development in computational 
processing power, agglomeration in fluid beds cannot be fully understood or modelled just by treating the 
binary coalescence situation, as it has been emphasised in Hede (2005 & 2006b). More advanced algorithms 
must be developed and implemented into soft-particle programs in order to account for the random nature 
of particle collisions in real fluid bed. Soft-particle simulation is nevertheless a promising tool for studying 
the effect at particle level of changing some of the physical parameters. Once the program is set up, changes 
can be made infinitely as only the processing time sets the limit with the present available computers. E.g. 
have the previously mentioned simulations by Lian et al. (1998) indicated that during binary collision the 
dominant energy dissipation is the viscous dissipation, except when the fluid viscosity is relatively small. 
These tendencies would have been extremely time consuming to extract from experimental data. 
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1.2.3 Summing up on Lagrangian modelling 

Since Lagrangian models describe particle motions in detail, it is expected that these models show closer 
resemblance with experimental results than with the Eulerian models. However, a direct comparison 
between hard- or soft-sphere models and experiments has not been made so far mainly because of the large 
number of particles that is required to justify the application of the continuum approach on one hand, and 
the limited number of particles that can be handled by the discrete element models on the other hand. 
Further complications arise from the fact that a rigorous comparison can only be made if the discrete 
element models account for the full three-dimensional motion of the particles as a two-dimensional 
modelling of the particle collision dynamics has proven to be too restrictive (Hoomans, 1999). This further 
strongly increases the required number of particles and consequently the computational demands. The 
number of particles that can be accounted for in such models is a generic but serious limiting factor in 
any of the present Lagrangian models. Even with modern computers, present models cannot account for 
more than 10 6 particles, which is several orders of magnitudes lower than that encountered in industrial 
fluid beds. Rough estimations based on Moores law for computer processing development estimate that 
models accounting for all three dimensional particles in industrial scale fluid beds will not be within 
reach in the coming ten to fifteen years (Webopedia, 2005 and Michales, 2003). 

1.3 Other granular dynamic modelling principles 

Besides the Eulerian and Lagrangian approaches, other principles have been applied for particle systems 
similar to fluid beds. These principles will only be briefly introduced, as the reported simulations of fluid 
bed systems using these techniques are extremely limited in number. 

1.3.1 Monte Carlo techniques 

Another method to study many particle systems is the Monte Carlo 7 technique. A Monte Carlo simulation 
is a mathematical experiment in which the behaviour of a system is simulated incorporating stochastic 
behaviour modelled using a randomness generator to vary the behaviour of the system (Kaye, 1997). 
This principle has been applied to several particle technology disciplines as reviewed by Wauters (2001). 
Regarding fluid bed granulation, Rosato et al. (1996) have applied Monte Carlo techniques. In their 
simulations a new overlap-free particle configuration was generated in each processing step. The change 
in the system energy was then calculated and if the change was negative, the new configuration was 
accepted. If the system energy on the other hand increased, the new configuration was accepted with a 
probability obtained from a statistical distribution based on the change in energy. With this method it 
was possible to simulate segregation and collision phenomena in agitated systems as fluid beds. 
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The Monte Carlo technique is capable of predicting steady state conditions being e.g. equilibrium 
conditions and for this purpose it has certain advantages over other hydrodynamic simulations 
techniques. As time is not a variable in Monte Carlo simulations, a Monte Carlo step can only be linked 
to an actual time step by means of calibration, which is a difficult task. In practice this means that Monte 
Carlo simulations of particle dynamics is not possible without input of a-priory knowledge. This is an 
obvious disadvantage when simulating agglomerating systems as the likelihood of permanent coalescence 
is closely related to time (Hoomans, 1999). Recent fluid bed hydrodynamic simulations seldom use 
the Monte Carlo principles although the technique has been applied elsewhere in fluid bed modelling. 
E.g. have Hapgood et al. (2004) used Monte Carlo techniques for the simulation of atomised droplets 
impacting the fluidised particle bed. 

1.3.2 Finite Element scheme 

The finite element method is well known from simple transport problems but has recently been applied 
to fluid bed systems. The principle involves discretising a large domain into a large number of small 
elements (which is often chosen to be triangles), developing element equations, assembling the element 
equations for the whole domain and then solve the assembled equations. The finite element discretisation 
of the governing differential equations is based on the use of interpolating polynomials to describe the 
variation of a field variable within an element. This makes finite element method well suited for irregular 
geometries and heterogeneous materials. Recent advances in finite element programs makes it possible 
to perform simulations that account for dynamic transport, heat and mass transfer, axial and radial 
dispersion, temperature and pressure variations and different hydrodynamic flow regimes, which is a 
prerequisite for fluid bed simulations (Wang & Sun, 2003). Results from catalytic fluid bed systems by 
Mahecha-Botero et al. (2005) seems promising as the models requires fewer assumptions than any of 
the principles presented earlier in this chapter. Modern finite element models can be reasonable accurate 
predictors of particle stresses if the bulk properties used as constitutive parameters are measured, but 
the finite element method struggles in general with boundary conditions when the bed particle material 
is fluidised at high fluidisation velocities (Bell, 2005). As it is the case with other simulation techniques, 
finite element model requires large processing times. 
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Summary 

The hydrodynamic approach on fluid bed granulation is a rapid growing field and a number of the most 
used available techniques for modelling the fluid bed particle-gas flow system have been presented in 
the previous sections. As in other parts of the chemical engineering science, there is a tendency to try 
to simulate and to model the hydrodynamics rather than to build uniform test equipment or try to 
conduct detailed experiments. With the complex behaviour of fluidisation and the chaotic phenomena of 
particles trajectories, the reasons for this are obvious. As the computational processing capacity is roughly 
doubled each 18 month, the boundaries for what is possible and what is not, are in continuous motion. 

Effective modelling of solid-fluid flow requires methods for adequately characterising the discrete nature 
of the solid phase and representing the interaction between solids and fluids. CFD multiphase models 
such as the Euler method address the problem within a continuum framework. In continuum models, 
contact between gas, particles and boundary surfaces is not considered explicitly with respect to particle 
inertia and mechanical properties. This limits the ability of CFD multiphase models to adequately 
represent particle-particle and fluid-particle interactions, and may therefore reduce the accuracy of the 
prediction of both the fluid and the particle dynamics. Despite the modelling challenges, applications 
of CFD to model hydrodynamics continues to develop as it has many advantages including design 
optimisation and scale-up of systems. Some of the correlations used in the present models however 
remain to be empirical or semi-empirical. As a result, the model and its parameters must be validated 
against experimental data obtained at similar scale and process configurations. The limitation of CFD 
models to represent particle-level interactions can be overcome by explicit calculation of the particle 
contact mechanics in a particle-scale reference frame using a Discrete Element (Lagrangian) approach 
such as the hard-sphere or soft-sphere principle. 

As it was presented earlier, a discrete element algorithm is basically a numerical technique, which 
solves engineering problems that are modelled as a large system of distinct interacting general shaped 
(deformable or rigid) bodies or particles that are subject to gross motion. Engineering problems 
that exhibit such large scale discontinuous behaviour as a particle fluid bed cannot be solved with a 
conventional continuum based procedure such as the Finite Element Method, although new expansions 
of the Finite Element technique in fact seems to have somewhat solved these issues. The discrete element 
procedure is used to determine the dynamic contact topology of the bodies. It accounts for complex 
non-linear interaction phenomena between bodies and numerically solves the equations of motion. 
Since the DEM is a very computationally intensive procedure, many existing computer codes are limited 
to model either two-dimensional or small three-dimensional problems that employ simple particle 
geometries. A general problem with Discrete Element modelling is that the simulation either concerns 
mechanical properties such as breakage and attrition or concerns agglomeration, but never accounts for 
all phenomena in the same model. 
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Despite the different theories and techniques associated with the presented hydrodynamics modelling 
principles, validation is a general problem with any of the simulation techniques. There are only few 
research groups around the world that work specifically in the field of fluid bed hydrodynamic modelling 
and most of the equipment, that is used to validate the simulation results, have been built by the research 
group itself. This means that published results are difficult to reproduce and further exploit and expand. 
As hydrodynamic modelling generally is very time-consuming and further requires advanced equipment 
such as high-speed cameras and specially designed fluid beds allowing the fluidisation behaviour to be 
recorded, hydrodynamic modelling will probably for some time on continue to be a part of fluid bed 
modelling only for an exclusive number of research groups around the world. As a sum-up, table 1 presents 
some of the most important advantages and disadvantages associated with the different hydrodynamic 
modelling techniques presented in the present chapter. 
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the coding part 

materials. 


demanding as the 

of particle collisions 

effects of particle 

of the simulation 



other techniques. 

and forces acting 

collisions and 

fast and reliable. 

FE models 



on the particles by 

forces acting on the 


require fewer 


CFD provides 

the gas. 

particles by the gas. 

The principles 

assumptions 


a detailed 



behind the 

than any of the 


understanding 

A larger number 

Soft-particle 

Monte Carlo 

other principles. 


of flow velocity 

of particles can be 

simulation is a 

techniques 



distribution, weight 

included into the hard 

promising tool for 

are easy to 

A range of 


loss, mass and 

sphere models than is 

studying the effect 

follow, and 

commercial FE 


heat transfer. 

possible in soft-sphere 

at particle level of 

the principles 

programmes 



models. 

changing some 

may be used 

is available, 


CFD makes it 


of the physical 

elsewhere in 

making the 


possible to fast 


parameters involved 

the granulation 

coding part of 


evaluate geometric 


in the granulation 

process - e.g. in 

the simulation 

(/) 

ai 

CT 

equipment 


process. 

the prediction 

fast and reliable. 

(0 

+■* 

changes with 



of the wetting 


<u 

> 

much less time and 


Theoretical particle 

behaviour of 


■o 

< 

cost than would 


level models may 

the atomised 



be involved in 


be validated using 

droplets on the 



laboratory or pilot 


the soft-particle 

fluidised particle 



plant testing. 


approach as 
numerous variations 

bed. 



A whole range of 


in the physical/ 

Monte Carlo 



CFD commercial 


chemical parameters 

models are 



programs are 


may be simulated 

capable of 



available making 


relatively fast once 

predicting steady 



the simulation 


the simulation 

state conditions. 



coding part fairly 
straightforward. 


program is set up. 






Soft-sphere models 




CFD can be applied 


are well suited 




in the process of 


for studying the 




scaling equipment 


modelling of 




as the CFD models 


impact breakage 




are based on 


of pre-existing 




fundamental 
physics and are 
thereby scale 
independent. 


agglomerates. 
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Summary 



Knowledge of the 

Hard-sphere models 

Soft-particle 

Limited available 

Limited available 


equation of state 

lack precision in fluid 

simulations 

literature makes 

literature makes 


for the particles is 

bed processes without 

struggle with high 

Monte Carlo 

Finite Element 


needed a-priori. 

rapid particle flow. 

computational 

simulations yet 

simulations yet 




processing 

an unproven 

an unproven 


Results from 

Although hard- 

demands. 

technique. 

technique. 


empty vessel CFD 

sphere simulations 





simulations may 

can yield both size 

Detailed information 

A-priori 

Problems with 


not readily be used 

and composition 

of binary collisions 

knowledge from 

boundary 


when particles are 

distribution of 

is far from being 

experiments is 

conditions 


being processed. 

granules, it is generally 

representative of the 

needed as input 

means that FE 



not suitable for 

situation inside a 

in the models. 

until recently has 


In CFD models 

realistic representation 

fluid bed. 


been unsuited 


for fluid bed 

of the granule 


As time is an 

for simulations 


granulation. 

microstructure (i.e. the 

Although soft- 

important 

of systems in 


contacts between 

internal distribution 

sphere simulations 

parameter in 

random motion. 


fluid, particles and 

of primary solids. 

can yield both size 

the granulation 


V) 

Qi 

boundary surfaces 

binder and porosity 

and composition 

process, time 

Large processing 

u\ 

(0 

are not considered 

of the granule). 

distribution of 

scale statements 

time even 

+■* 

c 

(U 

explicitly with 

Microstructure is an 

granules, it is 

should be 

with modern 

> 

-o 

respect to particle 

important property 

generally not 

incorporated 

computers. 

(0 

inertia and the 

especially in the case 

suitable for realistic 

into Monte Carlo 


5 

mechanical 

with enzyme granules 

representation 

models before 



properties of the 

as it determines the 

of the granule 

any advanced 



particles. 

release rate of the 

microstructure 

simulation of 




enzyme ingredient 

(i.e. the internal 

the fluid bed 



Due to the 

from the granule. 

distribution of 

granulation 



obvious need 


primary solids. 

process should 



for accounting 

Present models 

binder and porosity 

be reliable. 



precise particle 

are only capable 

of the granule). 




level properties 

of accounting for 

Microstructure 




into fluid bed 

50,000 granules at 

is an important 




hydrodynamic 

the time, thereby 

property especially 




models, CFD 

making simulations 

in the case with 




models must 

comparable only to 

enzyme granules as 




be regarded as 

experimental data 

it determines the 




inappropriate for 

from laboratory scale 

release rate of the 




simulating particle 

fluid beds. 

enzyme ingredient 




systems. 


from the granule. 




Table 1 : Advantages and disadvantages associated with the hydrodynamic modelling techniques presented in the text. 
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Table of symbols 


Nomenclature 


Unit (Si-system) 


a 

a' 

aa 

sae 

a n 

A 

A* 

b 

bAE 

C 

bair distrib pi. 

4 

dbed 

dd 

dd,rel 

dp 

dorifices 

dsp/sp 

d v 

dvessel 

e 

E 

E* 

Eelu 

fi(x, r, t) 
fbi 

^initial 

ftetra 

Fpend.,bound. 
Fpend.,eq sph. 
Fpend.,gorge. 
Fvis 


Fi 

Ft 

Fi E 

F E 


H 


Internal coordinate 
Material constant 

Projected area of liquid binder droplets 
Fitting parameter 

Projected area of a nucleus granule 
Powder flux 

Contact area between colliding granules 

Internal coordinate 

Fitting parameter 

Cohesivity of dry particle mass 

Air distribution plate diameter 

Gas bubble diameter 

Fluidised bed diameter 

Liquid droplet diameter 

Relative liquid droplet diameter 

Particle diameter 

Pitch orifice diameter 

Interaction parameter of two spheres 

Equivalent diameter of particles 

Fluid bed vessel diameter 

Particle coefficient of restitution 

Young modulus 

Granule Young modulus 

Elutriation rate 

Average number density function 
Bi-variant average number density function 
Initial average number density function 
Tetra-variant average number density function 
Pendular force in the “boundary” method 
Pendular force between two equally sized spheres 
Pendular force in the “gorge” method 
Viscous force 

Net force vector acting on particle i 
Drag force vector 

Force vector accounting for external fields 

Force vector accounting for particle-particle interactions 

Gravity 


Dimensionless 

m 2 

Dimensionless 

m 2 

m 2 /s 

m 2 


Dimensionless 

N/m 2 

m 

m 

m 

m 

m 

m 

m 

m 

m 

m 

Dimensionless 

N/m 2 

N/m 2 


N 

N 

N 

N 


m/s 2 
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G(x, r, t) 

Rate of growth by layering 

- 

G s 

Mass flux of particles 

m 2 /s 

h(x, r, t) 

Net generation rate of particles 

- 

ho 

Binder layer thickness covering colliding granules 

m 

ha + 

Birth of particles due to aggregation 

- 

h a ' 

Death of particles due to aggregation 

- 

hasp 

Characteristic length scales of surface asperities 

m 

h b + 

Birth of particles due to breakage 

- 

h b ' 

Death of particles due to breakage 

- 

hbed 

Bed height 

m 

H 

Separation distance between two spheres 

m 

i 

Summation parameter 

- 

h 

Moment of inertia 

- 

J 

Nucleation ratio 

Dimensionless 

k 

Proportionality constant 

- 

k' 

Proportionality constant 

- 

ken 

Coordination number 

Dimensionless 

K a 

Nucleation area ratio 

Dimensionless 

L 

Characteristic length of particles 

m 

Lbed 

Fluid bed length from distributor plate to exhaust exit 

m 

Lslr 

Length scale ratio 

Dimensionless 

m 

Mass 

kg 

m ( x ) 

Mass of a particle of internal state x 


m aggl 

Agglomerate mass 

kg 

IWbed 

Bed load 

kg 

I^harm 

Harmonic mean granule mass 

kg 

■^nozzle air 

Spray rate from the nozzle 

kg/s 

^ spray 

Air flow rate through the nozzle 

kg/s 

Mi 

Net torque vector 

- 

n fl 

Flow index 

Dimensionless 

n(x, r, t) 

Actual number density 

- 

h 0 

Nucleation rate 

No. of particles/s 

N(r, t) 

Average total number of particles 

- 

N t 

Total number of particles 

- 

Nto 

Initial total number of particles 

- 

P 

Summation number 

Dimensionless 

P 

Pressure 

Pa 

P(x,r|x',r) 

Probability density function 

- 

q 

Discretisation number 

Dimensionless 

r 

Radius 

m 

r aggl 

Radius of an agglomerate 

m 
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* 

rdef 

Critical radius of an agglomerate after which def. occurs 

m 

^harm 

Harmonic mean granule radius 

m 

Tneck 

Pendular bridge neck radius 

m 

F g 

Mean granule size 

m 

F g» 

Initial mean granule size 

m 

Tieck 

Pendular bridge neck radius 

m 

r 

External coordinate vector 

- 

r' 

External coordinate vector 

- 

R 

Radius 

m 

R P 

Particle radius 

m 

®max 

Maximum pore liquid saturation 

Dimensionless 

s 

Distance 

m 

Sc 

Saturation at transition funicular/capillary state 

Dimensionless 

Sd 

Dry coating material feed rate 

- 

Sf 

Saturation at transition pendular/funicular state 

Dimensionless 

Ssat 

Amount of saturation 

Dimensionless 

Stdef 

Stokes deformation number 

Dimensionless 

st de ; 

Critical Stokes deformation number 

Dimensionless 

St v 

Viscous Stokes number 

Dimensionless 

St v * 

Critical viscous Stokes number 

Dimensionless 

Sw 

Wetting saturation 

Dimensionless 

S(q) 

Summation function 

- 

S Kolmogorov 

Kolmogorov entropy 

bits/s 

t 

Time 

s 

fcoat 

Coating time 

s 

U 

Granule velocity 

m/s 

U 0 

Initial granule collision velocity 

m/s 

u 

Fluidisation velocity 

m/s 

u br 

Bubble rise velocity for a fluid bed 

m/s 

u mf 

Minimum fluidisation velocity 

m/s 

U s 

Superficial gas velocity 

m/s 

Vi 

Velocity vector 

- 

V 

Particle volume internal coordinate 

- 

V 

Average particle volume 

m 3 

Vl 

Liquid binder volume internal coordinate 

- 

V 

Volumetric spray rate 

m 3 /s 

Vaggl 

Agglomerate volume 

m 3 

V bridge 

Liquid bridge volume 

m 3 

v r 

Volume of external coordinates 

- 

V x 

Volume of internal coordinates 

- 

w 

Granule volume parameter in coal, kernel expression 

- 
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* 

w 

Critical average granule volume 

- 

W m 

Mass ratio of liquid to solid 

Dimensionless 

W 

Spray zone width 

m 

X 

Internal coordinate vector 

- 

x' 

Internal coordinate vector 

- 

X 

Coordinate 

m 

y 

Coordinate along the width of the spray zone 

m 

Y(r,t) 

Continuous phase vector 

- 

Y d 

Plastic yield stress 

N/m 2 

z 

Counting number 

Dimensionless 
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Greek 


p 

Coalescence kernel 

- 

Po 

Rate constant 

- 

Pdt 

Aggregation probability in time interval dt 

- 

Pid 

Coefficient of interphase drag 

Dimensionless 

P‘ 

Coalescence kernel expression 

- 

5 

Dimensionless bubble spacing 

Dimensionless 

bpdef 

Extent of permanent plastic deformation 

Dimensionless 

Ff 

Coefficient of internal friction 

Dimensionless 

Fmean 

Mean in the Gaussian distribution 

m 

Of 

Macroscopic shear stress at failure 

Pa 


Macroscopic normal stress 

Pa 

^width 

Standard deviation 

m 

G t,f 

Funicular bridge static tensile strength 

N/m 2 

^t,p 

Pendular bridge static tensile strength 

N/m 2 

°t,c 

Capillary bridge static tensile strength 

N/m 2 

c(y) 

Characteristic stress in an agglomerate 

N/m 2 

a y 

Yield stress/strength 

N/m 2 

X c 

Average particle circulation time 

s 

Xd 

Droplet penetration time 

s 

Va 

Dimensionless spray flux 

Dimensionless 

Vn(y) 

Dimensionless nuclei distribution function 

Dimensionless 

Vn 

Dimensionless spray number 

Dimensionless 

C 

Particle shape factor (sphericity) 

Dimensionless 

l 

Dimensionless parameter in the dynamic strength eq. 

Dimensionless 

<P 

Half filling radius 

o 

0 

Contact angle 

o 

V 

Poisson ratio 

Dimensionless 

v(x',r',Y,t) 

Average number of particles formed from break up 

- 

8 

Particle voidage (void fraction) 

% 

^longitudinal 

Longitudinal extension strain 

Dimensionless 

£min 

Minimum porosity 

% 

^trans 

Transverse contraction strain 

Dimensionless 


Mean granule porosity (void fraction) 

% 

Yiv 

Interfacial surface tension between liquid and vapour 

N/m 

Y 

Shear rate 

s" 1 

P 

Density 

kg/m 3 

Pb 

Binder liquid density 

kg/m 3 
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Pg 

Granule density 

kg/m 3 

Pp 

Particle density 

kg/m 3 

'lapp 

Apparent viscosity 

kg s /m 

llliq 

Liquid (binder/coating) viscosity 

kg s /m 

COi 

Angular velocity vector 

- 

Q h 

Hounslow discretisation parameter 

- 

Q r 

Domain of external coordinates 

Dimensionless 

a 

Domain of internal coordinates 

Dimensionless 
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1. Some of the newest studies in granulation modelling do however indirectly include population balance 
equations as well as discrete element simulation techniques (e.g. Gantt & Gatzke, 2005). 

2. Fluid dynamics is the sub discipline of fluid mechanics that studies fluids in motion. Fluids are specifically 
gases and liquids. Commonly fluid dynamics is divided into sub disciplines being aerosol dynamics (the 
study of gases) and hydrodynamics (the study of fluids and thereby fluid like systems) (Wikipedia, 2005). 

3. Also known as “Continuum models” or “Eulerian-Eulerian models” (Goldschmidt, 2001 and Taghipour 
et al., 2005). 

4. Actually Computational Fluid Dynamics is formally the overall term for all hydrodynamic modelling 
(and thereby also for DEM), but it has become common to refer to other terms when granular dynamics 
is involved (Hoomans, 1999). 

5. An excellent review of available commercial CFD programmes can be found in Xia & Sun (2002). 

6. Hoomans (1999) has compared the results of a soft-sphere model to those of a hard-sphere model and 
found that, provided that the spring stiffness was high enough, the difference in bed hydrodynamics 
between the two types of models were very small. This further supports the soft-particle approach for 
fluid bed modelling. 

7. The term “Monte Carlo” refers in general to any technique of statistical sampling employed to approximate 
solutions to quantitative problems such as dice tosses or the chance that two particles in a bulk sample 
will collide (Webopedia, 2005). 
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